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Using exact diagonalizations, we investigate the T = 
phase diagram of the Multi-Spin Exchange (MSE) model on 
the triangular lattice: we find a transition separating a fer- 
romagnetic phase from a non-magnetic gapped Spin Liquid 
phase. Systems far enough from the ferromagnetic transition 
have a metamagnetic behavior with magnetization plateaus at 
m/m 3 at = and 1/2. The MSE has been proposed to describe 
solid 3 He films adsorbed onto graphite, thus we compute the 
MSE heat capacity for parameters in the low density range of 
the 2 nd layer and find a double-peak structure. 

PACS numbers: 75.10.Jm; 75.50.Ee; 75.40.-s; 75.70.Ak 
An increasing number of experiments on 3 He films have 
enforced the triangular lattice Multiple-Spin Exchange 
(MSE) picture of the solid second layer ( [|l]j^] and ref- 
erences therein). Following Thouless the magnetic 
properties of a 2-dimensional quantum crystal are de- 
scribed in spin space by the effective Hamiltonian: 
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where P is any permutation operator of the spins of the 
lattice, and Jp equals one half of the (positive) tunneling 
frequency associated to the exchange process P. Whereas 
exchange of an even number of spins favors antiferromag- 
netism, odd processes are ferromagnetic. On the trian- 
gular lattice with spin 1/2, two- and 3-body exchange 
reduce to an Heisenberg Hamiltonian with an effective 
= Ji — 2 J3. In two-dimensional 3 He, because of sto- 
ichiometric hindrance, the 3-body process is much more 
efficient than the 2-body one, and the effective coupling 
constant is ferromagnetic. However, it was suspected 
since a long time that n-particle terms with n > 3 could 
not be ignored 0,0,0] ■ Very recently, a thorough analysis 
of susceptibility and heat capacity measurements enabled 
Roger and the Grenoble group (7) to establish the density 
dependence of the J's and the importance of the 4-spin 
exchange J4 in a large density range of the 2 nd layer: from 
the 2 nd layer solidification (AF solid) to the third layer 
promotion (F solid). This is in accordance with recent 
Path Integral Monte Carlo calculations which estimate 
Jf/Ji ~ -2± 1, J 5 /J 4 ~ 1 ± .5 and J 6 /J 4 ~ 1 ± .5 §. 
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As first suggested by Roger in 1990 |J, 4-spin exchange 
strongly frustrates the system. Momoi, Kubo and Niki 
have recently studied the J| ff -J4 model in the classical 
and semi-classical (spin- wave) limits and found numerous 
ordered phases: a ferromagnetic, a 4-sublattice ferrimag- 
netic phase, 3- and 4-sublattice AF phases and a chiral 
ordered phase Our SU(2) Schwinger-Boson anal- 

ysis of the Jf ff - J4- J5 Hamiltonian also pointed to a rich 
phase diagram with Neel as well as helicoidal phases [[l2| . 

In the present work, exact diagonalizations results 
show that none of these T=0 Long Range Ordered (LRO) 
AF phases survive spin-1/2 quantum fluctuations. In- 
stead, in the AF region, the system is a Quantum Spin 
Liquid (QSL) with short range spin-spin correlations, as 
suggested by Ishida et al. |Q and Kubo et al. jll} . 

We truncate Eq. ([j]) to the following simplest cyclic 
exchange patterns: 

H=jfY J P* + J* £ + (2) 
-J 5 £ (P 5 + P 5 - 1 ) + J e £ (Pa + Pf 1 ) 

^ O 

The spin-1/2 permutation operators can be rewritten 
with usual spin operators (Pauli matrices): Py = 
2Sj.Sj + 1/2. P1234 + h.c. and P12345 + h.c. (resp. 
Pi23456 + ^-c) are polynoms of degree two (resp. three) in 
Si.Sj [|l2|jl4| . The quantum phase diagram of this model 
is studied through the analysis of the finite size scaling 
of the low energy spectra of samples subjected to various 
boundary conditions (twisted or periodic with different 
shapes). 

Ferro-antiferro transition. — We first look at the T = 
ferro-antiferro transition line of Eq. (0). This line (see 
Fig. |]) is determined from about 50 spectra (N = 19) 
in the J n space, with most points near the transition. 
Far all cases, the ground-state is either an S = or 
an S — N/2 state. This excludes the possibility of a 
uuud phase found in the classical calculations of Kubo 
et al. [ [TT| |, which is ferrimagnetic and has a total spin 
S = N/4. The line where J x , the 1/T 2 coefficient of 
the susceptibility, vanishes stands roughly parallel to the 
T = F/AF line, inside the AF region. The density 
dependence of the J n 's proposed by Roger et al. for the 
second layer (see crosses in Figure |l|) leads us to conclude 
that a T=0 transition to ferromagnetism occurs at P2 — 
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6.8 ± 0.3nm -2 whereas J% is zero at pi ~ 6.5nm~ 2 0. 

iVo Long Range Order. — The nature of the non mag- 
netic or antiferromagnetic (AF) phase is a more challeng- 
ing question. We first look for signatures of Neel Long 
Range Order (NLRO). NLRO is characterized on finite 
systems by very stringent spectral properties |15|: i) The 
symmetry breakings associated to the order parameter 
are embodied in a family of ~ N a low lying levels col- 
lapsing to the ground-state in the thermodynamic limit 
(a is the number of magnetic sublattices). These levels 
with definite space and SU (2) symmetries and dynamical 
properties should appear directly below the first magnon 
excitations, ii) The finite-size scalings of these levels are 
known. In particular, the ground-state energy per site 
E(S — 0,N)/N has corrections scaling as iV~ 3 / 2 , the 
AS = 1 spin gap E(S — l,N) — E(S = 0,N) goes to 
zero as N" 1 . None of these prescriptions is obeyed by 
the MSE spectra in the "S=0 ground-state" region dis- 
played Fig. [I], whatever the twisted or shape boundary 
conditions may be. Thus, we exclude any commensurate 
or non-commensurate NLRO. 

A — Ji model. — We have analyzed in detail the 
= — 2, J4 = 1 model for nearly all possible systems 
from N — 6 to 30 and for N — 36. We noticed two 
different scaling behaviors: samples with N multiple of 
4 or 6 have a low ground-state energy increasing with N 
whereas others samples have a high ground-state energy 
decreasing with N. The energies of both families merge 
for No ~ 40. Our interpretation is the following: No is 
a crossover size above which the system is not anymore 
sensitive to boundary conditions and this is the signature 
of a finite length scale £ = \/Nq in the ground-state wave 
function. This is supported by two facts: a fast decay 
of the spin-spin correlations with distance and a non- 
vanishing spin gap A in the thermodynamic limit. AS = 
1 gaps are plotted in Figure |J the two families of samples 
(squares for N multiple of 4 or 6 and stars for others) 
have gaps of the order of 1 for N > 24. An estimation of 
the N = 00 gap is possible using the strong correlation 
between E/N and A, the result is A ~ 1.1 ± 0.5 (details 
will be given elsewhere). 

These data point to a Quantum Spin Liquid state with 
a gap of the order of 1 for J| ff = — 2 and J4 = 1. This gap 
and the sensitivity to the geometrical shapes and bound- 
ary conditions of the small samples suggest a valence- 
bond picture of the ground-state and of the first triplet 
excitation. 

Because of the strong frustration between the effec- 
tive first neighbor ferromagnetic Heisenberg term (Jf ) 
and the 4-spin antiferromagnetic exchange, it is the tri- 
angular 6-spin plaquette which is the first system with 
a paramagnetic (S — 0) ground-state and a significant 



gap[]. However the ground-state wave function is not a 
naive tensorial product of S = independent triangle 
wave functions: such an approximation gives a very high 
ground-state energy (—2.8 to be compared to —4) and 
largely underestimates the gap (.4 to be compared to 1). 
This quasi independent triangle picture is thus deeply 
renormalized by resonances. 

Analysis of the low lying levels in the S — subspace 
leads to the same conclusion: the number of singlets 
below the first triplet level is very small (< 10 all de- 
generacy taken into account). As a comparison, this 
is very different from the Kagome case where a con- 
tinuum of singlets state is found in the magnetic gap 
]l6|,0. Thus, this system seems a very example of a 
short range Resonating Valence Bond state with a clear 
cut gap in a translationally invariant 2-dimcnsional spin- 
1/2 model. Consequently, both the low temperature spe- 
cific heat and spin susceptibility arc thermally activated 
(x(T) ~ g— a/t^ c v (T) ~ e' A / T ). However the gap de- 
creases rapidly when approaching the AF/F transition 
and we are not yet able to decide if the gap vanishes 
at the transition to ferromagnetism or before. Indeed, 
experimental results of Ishida and collaborators seem to 
indicate either a very small gap or a quantum critical 
behavior |13|. 

,2 nd layer heat capacity. — We compute the MSE heat 
capacity for J 2 cff /J 4 = -2, J5/J4 = 0.2, J 6 /J 4 = 0.08 
with N = 16,20 and 24 (Fig. |) and find a clear 
low temperature peak. The entropy at T/Jc v — -5 is 
~ 0.4JV7n(2) and the low temperature peak is thus likely 
to remain in the thermodynamic limit. It also subsists 
for a relatively large range of competing J n . Such a low 
temperature peak is characteristic of different highly frus- 
trated systems |l^||- The high temperature peak is lo- 
cated at a temperature of the order of Jc v (Jc v is the 
leading coefficient of the 1/T expansion of the specific 
heat: C v = 9/4{J Cv /T) 2 + OCl/T 3 ), its expression as a 
function of the J's is given in [19[). The low temperature 
peak height and location do not only depend on Jc v but 
also on the relative values of the coupling parameters. A 
better agreement between present results and experimen- 
tal results jl3| is expected by a fine tuning of the coupling 
parameters. Indeed, moving towards the boundary line 
between AF and F phases both decreases the gap and 
shifts the low temperature peak towards lower tempera- 
ture. 

Low energy degrees of freedom. — To understand the 



For N = 24 (resp. N = 12) one can compare the spec- 
trum of the sample built with four (resp. two) 6-site triangles 
with the spectra computed from other shapes. The triangle- 
compatible shape gives the lowest energy and the largest gap 
(others have energies and gaps comparable with the frustrated 
family ones). 
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excitations responsible for this low energy peak we looked 
for possible common properties of low lying levels: it ap- 
pears that most of these levels have a significant pro- 
jection on the subspace engendered by spin-1 diamond 
tilings. This subspace £ is defined by the non-orthogonal 
family of wave functions: 

l*>= (g) \S d = 1, SI G {-1,0,1}) (3) 

where \S d = l,S%) is the S = 1, S z — S% state, symmet- 
ric with respect to the small diagonal of the d th diamond. 
For N = 16 the projections of the exact low lying levels 
on £ range from 10% to 50% for nearly all the states in- 
volved in the peak (and drop under 1% for higher energy 
states) . These numbers are very large compared with the 
expectation values of the projection of a random 5 = 
or S = 4 wave function, which are, for the same lattice 
size, of the order of 1%. Since each tiling of the lattice 
gives 3t independent states, the entropy associated with 
£ is a least in(3"* ) ~ 0ANln(2), in agreement with the 
low temperature peak entropy found in our samples. 

These results lead to the following picture: At high 
temperature down to T ~ Jc v the degrees of freedom 
are essentially random spin- 1/2. For T less than Jc v , the 
thermal wavelength increases and near neighbor spins be- 
have coherently as weakly-ferromagnetic entities (pseudo 
spin-1). This explains the low temperature peak. At 
T ~ A (spin gap), the 4-spin exchange coupling creates 
larger clusters and the system organizes itself as a QSL. 

Magnetization. — Now we apply a magnetic field and 
look for the spin S of ground-state (T=0 magnetiza- 
tion). The large gap between the sectors of total spins 
S = S max /2 and S = S max /2 + 1 indicates that excit- 
ing one diamond to its S = 2 state costs a large energy. 
This feature gives rise to a low temperature plateau^] at 
magnetization m = 1/2 (Fig. ||) which has also been 
found in the classical variational picture by Kubo et al. 
p| . For the same coupling parameters as Fig. ||, the 
figure ^ shows two metamagnetic transitions at Hex and 
Hc2 and the magnetization is completely quantized^] : 
to = 0, 1/2 or 1. Thanks to the large gaps (S = — > 1 
and S = N/4 — » N/A + 1) these plateaus survive thermal 
excitations. Close to the ferromagnetic transition (upper 



2 Related phenomena have already been encountered in the 
Heisenberg model with Ising anisotropy in the MSE 
model on the square lattice ]2lj], on Heisenberg ladders or 
chains |^,^3). An m = plateau has been observed in the 
spin- ladder compound Cu 2 {C 5 H w N 2 )2Ch 

3 This, as well as the symmetries and degeneracies of the 
low lying levels of the spectra, will be discussed elsewhere 
in relation with possible extension of the Lieb-Schultz-Mattis 
theorem to 2-dimensional magnets |26|27j]. 



black triangle of Fig. |lj) both gaps decrease and the over- 
all shape of the energy versus magnetization anticipates 
a transition to a phase where the magnetization curve is 
strongly non linear and reminiscent of the metamagnetic 
transitions (see [pl| about metamagnetism in 3D solid 
3 He). 

In this letter we have shown that the spin-1/2 MSE 
model, in the range of parameters relevant for the de- 
scription of 3 He on graphite, exhibits a transition from 
an antiferromagnetic phase at low density to a ferromag- 
netic one at high density. The AF phase has no Neel 
Long Range Order and is a Quantum Spin Liquid phase. 
In this phase the specific heat has two peaks. The low- 
est one testifies the building of a QSL out of short range 
pseudo spin-1. This picture is consistent with magnetiza- 
tion and heat capacity measurements. The microscopic 
arrangement of the spins may also show up as plateaus in 
the magnetization and T = metamagnetic transitions 
between zero, half polarized and fully polarized phases. 
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FIG. 1. Phase diagram for jf < 0. The solid line is the 
T — transition line between ferromagnetic and antiferro- 
magnetic phases. The crosses are Roger et aFs results [jjj for 
four 2 nd layer densities (nm -2 ): p% = 6.5 (Je/Jn — 0.7), 
p\ = 7.0, pi = 7.65 and p d 2 = 7.8. For b, c and d Ja/Ji ~ 0.4 
(the size of the crosses may underestimate the uncertainties 
in the parameters). Two of the black triangles indicate the 
sets of coupling parameters corresponding to Fig. The 
upper one is the point close to the frontier mentioned in the 
text. 
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Spin gap plotted as a function of 1/N for 
2, J 4 — 1. The vertical bar is the N = oo ex- 
trapolation made out of the E/N <-> A correlation analysis. 
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3. Heat capacity versus temperature for three 



sizes. Coupling parameters are Jf /J4 = —2, J5/J4 = 0.2, 
Jul J a = 0.08. For these values Jc„ = 0.93J4 and the N — 00 
spin gap is of the order of J4/2. Due to finite-size effects, on 
N — 20 and 24 the high temperature peak (T ~ Jc v ) only 
shows up as a shoulder. 
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FIG. 4. Magnetization versus magnetic field B. Exchange 
parameters are those of Fig. ^. Since in 3 He, J4 ~ 2mK and 
g/i = 1.5mK.T~ , Hci is in the Telsa range. 



4 



